Single cell studies of the cell cycle and some models

Analysis of growth and division often involves measurements made on cell populations, which tend to average data. The value of single cell analysis needs to be appreciated, and models based on findings from single cells should be taken into greater consideration in our understanding of the way in which cell size and division are co-ordinated. Examples are given of some single cell analyses in mammalian cells, yeast and other microorganisms. There is also a short discussion on how far the results are in accord with simple models.     

An approach to large scale identification of non-obvious structural similarities between proteins

Background

A new sequence independent bioinformatics approach allowing genome-wide search for proteins with similar three dimensional structures has been developed. By utilizing the numerical output of the sequence threading it establishes putative non-obvious structural similarities between proteins. When applied to the testing set of proteins with known three dimensional structures the developed approach was able to recognize structurally similar proteins with high accuracy.

Results

The method has been developed to identify pathogenic proteins with low sequence identity and high structural similarity to host analogues. Such protein structure relationships would be hypothesized to arise through convergent evolution or through ancient horizontal gene transfer events, now undetectable using current sequence alignment techniques. The pathogen proteins, which could mimic or interfere with host activities, would represent candidate virulence factors.The developed approach utilizes the numerical outputs from the sequence-structure threading. It identifies the potential structural similarity between a pair of proteins by correlating the threading scores of the corresponding two primary sequences against the library of the standard folds. This approach allowed up to 64% sensitivity and 99.9% specificity in distinguishing protein pairs with high structural similarity.

Conclusion

Preliminary results obtained by comparison of the genomes of Homo sapiens and several strains of Chlamydia trachomatis have demonstrated the potential usefulness of the method in the identification of bacterial proteins with known or potential roles in virulence.

     

A GWAS assessment of the contribution of genomic imprinting to the variation of body mass index in mice

Background

Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis.

Results

10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10⁻⁴ to 1.842 ×10⁻⁴ when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered.

Conclusions

When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).     

Defining window-boundaries for genomic analyses using smoothing spline techniques

Background

High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome.

Results

Simulations applying this method were performed to identify selection signatures from pooled sequencing F_ST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing F_ST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach.

Conclusions

We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on F_ST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.